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Summary. We discuss an effect of the electrostatic field on superconductivity near 
the surface. First, we use the microscopic theory of de Gennes to show that the 
electric field changes the boundary condition for the Ginzburg-Landau function. 
Second, the effect of the electric field is evaluated in the vicinity of -ff C 3, where the 
boundary condition plays a crucial role. We predict that the field effect on the surface 
superconductivity leads to a discontinuity of the magnetocapacitance. We estimate 
that the predicted discontinuity is accessible for nowadays experimental tools and 
materials. It is shown that the magnitude of this discontinuity can be used to predict 
the dependence of the critical temperature on the charge carrier density which can 
be tailored by doping. 



1.1 Introduction 

The surface of a superconductor is an important region in which the super- 
conductivity nucleates and which represents a natural barrier for penetrating 
or escaping vortices. It is desirable to control surface properties so that the 
nucleation can be stimulated or suppressed. An even more attractive task is 
to open or shut the penetration barrier for vortices. A promising tool of the 
surface control is the gate voltage for which we can benefit from the extensive 
technological experience with field effect transistors. 

Unfortunately, the interaction of the electric field applied to the metal 
surface with the superconducting condensate is very weak. Indeed, the su- 
perconducting condensate does not interact with the electrostatic potential 
as shown by Anderson [2]. The condensate feels only the indirect effects like 
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changes of the local density of states or eventual changes of the surface crystal 
structure. 

It is very likely that it will become possible to enhance the field effect on 
the superconductivity by a proper surface treatment. To this end it would be 
of great advantage to understand how the field interacts with the condensate 
and to have reliable experimental methods directly aiming to measure the 
strength of this interaction. 

To support the experimental effort in this direction, in this chapter we pro- 
vide a phenomenological theory of Ginzburg-Landau (GL) type supplemented 
with the de Gennes boundary condition derived from the microscopic Bardeen, 
Copper, and Schricffcr (BCS) theory. It will be shown that the boundary con- 
dition captures the field effect on the condensate while the GL equation deter- 
mines how the condensate responds to the field-affected boundary condition. 

The field effect on the superconductivity has been measured under various 
conditions, nevertheless its actual strength is not yet accurately established. 
The most pronounced field effects are observed on thin layers, in which it is 
possible to increase or lower their critical temperature [1,7,8,19,28]. These 
samples are so thin that the applied field considerably changes the total den- 
sity of electrons and the observed effect can be interpreted in terms of the 
modified bulk properties. With thicker samples one meets the problem that 
the potential field effect is restricted to the surface and the underlying bulk 
overrides its contribution. 

At the end of this chapter we discuss the field effect on the surface su- 
perconductivity near the third critical magnetic field H c ^. In this regime the 
bulk superconductivity is absent and the surface superconductivity crucially 
depends on the boundary condition. We will show that the field effect can be 
observed via the discontinuity in the magnetocapacitance [20,21]. 

1.2 Limit of large Thomas-Fermi screening length 

To introduce the field effect on the superconductivity we start from the theory 
employed by Shapiro and Burlachkov [3,23-27] and by Chen and Yang [4]. It 
is justified for high-T c superconductors in which the GL coherence length £ is 
very short, while the hole density is low leading to relatively large Thomas- 
Fermi screening length Atf- In these materials Atf ~ £, which allows us to 
introduce field induced effects via local changes of the parameters of the GL 
theory. 

Let us assume the jelly model in which the electric charge of electrons is 
compensated by a smooth positively charged background. Both charges are 
restricted to the half space x > 0. The electric field applied to the surface is 
exponentially screened E(x) = E e _x / ATF inside the metal. According to the 
Gauss equation edivi? = p, the induced electron density Sn — p/e reads 

Sn(x) = 4^e" x / ATP . (1.1) 
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In the GL equation 
1 



2m* 



(-iftV - e*Ay V + aip + P\ip\ 2 ip = 0, (1.2) 



the change of the electron density leads to changes of the GL parameters a 
and (3. According to the Gor'kov theory [9] (3 = 7 ^( 3 ) 1 | F <: Tt ■ It is thus a robust 
parameter in which Sn can create relative changes of the order of 8(3/(3 ~ 8n/n. 
These changes can be neglected. The other parameter a — a'(T — T c ) is 

a difference of two large constituents. Here, a' — 7 ^^ 3 f E ° is modified also 
only negligibly as 8a' /a' ~ 8n/n. But, for temperatures close to the critical 
temperature, T — > T c , even small changes in T c lead to large relative changes 
of a eventually even changing its sign. The GL equation with the dominant 
part of the field effect thus reads 

1 BT 

{-iUV - e*A) 2 4> + aip- a'^Snip + (3\ip\ 2 ^ = 0. (1.3) 



2m* dn 

We have denoted the field affected GL function as ip. Initially we shall use the 
customary GL condition J^tp = for this equation. 

Let us assume now that this theory holds also for conventional supercon- 
ductors where one always finds a sharp inequality Atf <C £■ Below we confirm 
the result obtained with this unjustified assumption by using the well justified 
microscopic approach of de Genncs. 

We split the GL function according to tp = ip + Sip, where Sip(x) = 
Sipe~ x ^ TF is the part of field induced perturbation which changes on the 
short scale Atf and tp covers the rest changing on the larger scale £. Our aim 
is to establish approximative <5?/> and to eliminate it, so that in the second step 
we will be left with the GL equation for the slowly varying function tp. 

The short scale component has an enormously large space gradient which 
dominates its contribution to the GL equation, 

: V 2 StP(x)-a'^5n{x)iP(x)^0. (1.4) 



2m* dn 

Since in this approximation the function 8tp is nonzero only in the narrow 
layer x <~ Atf , we can neglect the space dependence of ip and use its surface 
value. Performing derivatives one finds 

2m*A TF ,dT c eE , . . 

8ip = -^a'/— # 1.5 

h on e\TF 

Besides the very local contribution expressed by the function 8ip, the elec- 
tric field induces also a perturbation on the scale of the GL coherence length 
£. Indeed, the GL boundary condition demands the zero derivative of the total 
function 

|-(V, + <tyO = o (1.6) 
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and the local part has a nonzero derivative 

From the GL boundary condition (1.6) and relations (1.5) and (1.7) one thus 
finds the boundary condition for ip 



chp_ 
dx 



x=0 



2m* ,dT c eE , _ 

^--^(0). (1.8) 



This relation can be interpreted as a field-affected GL boundary condition. 

In usual problems handled by the GL theory one ignores the exact gap 
profile at the surface and focuses on its behavior deeper in the bulk on the 
scale £. In the same spirit we can ignore the short scale component Sip using 
the approximation ip w tp. Doing so, we have to keep in mind that the short 
scale component leads to the field-affected boundary condition (1.8). Solving 
the GL equation (1.2) with the boundary condition (1.8) one obtains the GL 
function where the field effect manifests itself on the scale £. 



1.3 de Gennes approach to the boundary condition 

The limit of large screening length does not apply to conventional supercon- 
ductors. In fact, in metals the Thomas-Fermi screening is smaller than the 
interatomic distance and the jelly model is not justified to describe the inter- 
action of the surface with the electric field. Naturally, the gradient correction 
represented by the kinetic: energy of GL theory is not a sufficient approxima- 
tion of the non-local part of the BCS interaction kernel. 

To access short screening lengths, Shapiro and Burlachkov [3,23-27] have 
employed a more sophisticated version of the GL theory in which the 'ki- 
netic energy' is not a mere parabolic function of the gradient but includes all 
derivatives up to infinite order in a form of the di-gamma function. Since high 
order gradients are important only in the short scale component, the above 
approach can be easily modified in this way. One merely replaces the kinetic 
energy in (1.4) by the corresponding di-gamma expression. 

Although this high order gradient correction is elegant and simple, it is 
likely not sufficient to cover changes on the sub- Angstrom scale, i.e., on the 
scale typical for majority of metals. Apparently, one should employ the mi- 
croscopic approach of Bogoliubov-de Gennes type pioneered by Koyama [13] 
and other groups [11,15-17,29-32]. Beside microscopic details of the gap near 
the surface, one has to take into account that the simple exponential decay 
of the charge from a sharp surface is not a very realistic model of metals, 
because electrons tunnel out of the metal. Realistic studies of the surface are 
problematic, however, even in the normal state, namely due to the nontrivial 
exchange-correlation interaction at the strong gradient of the electron density 
innate to all surfaces. 
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Instead of improving the above approach, it is advantageous to formulate 
the boundary condition directly from the microscopic theory of de Gennes [5] . 
De Gennes did not assume the electric field explicitly. His result, however, 
does not specify the forces forming the surface so that the applied field can 
be included. 

Gor'kov has shown that the BCS gap A and the GL function are propor- 
tional to each other, ip = const x A Using an extrapolation of the BCS gap 
from the bulk towards the surface, de Gennes has arrived at the boundary 
condition of the form 



1 dip 



V>(0) dx 



x=0 



1 OA 
A dx 



x=0 



(1.9) 



where b is called the extrapolation length. 

Within the BCS theory, the extrapolation length is given by the formula 



e(o)N v 



OO 

oV J 



dx 



A(x) 



1 - 



N(x) 



(1.10) 



derived by de Gennes (Eq. (7-62) in Ref. [5]). Here N(x) is the local density 
of states at the Fermi level and N is its bulk limit. The notation of the gap 
is different. While A(x) is the true local value of the BCS gap, A is the fake 
surface value after all short scale components have been removed, i.e., it is 
a value extrapolated from the near vicinity to the surface. Finally, V is the 
BCS interaction and £(0) is the GL coherence length at 'zero' temperature. 
In pure metals it is linked to the BCS coherence length £o as £(0) = 0.74£o- 

De Gennes estimated a typical value of b ~ 1 cm at metal surfaces in 
vacuum. This value is very large on the scale of the GL coherence length, 
therefore this contribution is usually neglected. The approximation, 1/6 pa 0, 
corresponds to the original GL condition J^ip = 0. 

Our aim is to include the effect of electric fields on the extrapolation 
length b. We denote as b the extrapolation length in the absence of the 
applied electric field and 5(l/b) = 1/6 — 1/bo reflects variation of the inverse 
length. The electrostatic potential corresponding to the electric field modifies 
the potential profile near the surface. It results in a change of the density of 
states SN(x). The density of states affects the gap function and creates its 
deviation 5A(x). In the linear approximation from (1.10) we find the change 
of the inverse extrapolation length as 



1 



em n v 



dx 



SAjx) 
A 



N(x) 
N 



A(x) 6N(x) 
A N 



(1.11) 



To estimate this change we recall the local density approximation in which 
the local density of states is a function of the local density, N(x) = N[p(x)]. 
Since the charge density 5p(x), which screens the applied electric field, spreads 
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over a layer of a few Angstroms near the surface, the perturbed density of 
state is restricted to this very narrow layer, too. We can thus neglect the 
^-dependence of A{x) in the second term and write 



){ b) e(0)N Q V 



dx 



5A(x) 



N{x) 



A(0) 6N& 
~A~o W 



(1.12) 



where 



5N^ 



OO 

■/ 



dx 6N(x) 



(1.13) 



is the total change of the density of states per area. 

To estimate the second term we assume that the local density of states 
achieves the bulk value on the scale of a few Angstroms. Since the gap function 
changes on the scale of the BCS coherence length which is much larger, we 
can expect that in the region of non-zero function 1 — N(x)/Nq the gap func- 
tion keeps its shape. Assuming that SA(x) w A(x)Sc, the last x-integration 
becomes identical to the integral in the de Gennes condition (1.10) so that 
(1.12) simplifies to 



5c 



1 



1 A(0) 5N (2 *> 



e(0)N o V A N 



(1.14) 



The relative change of the gap 5c can be estimated from the GL theory. 
The GL function obtained with the boundary condition of a large but finite 
extrapolation length b can be written as a sum of the constant bulk term 
V'oo = \f—uj~& and a small perturbation ip', i.e. ip = -0QO + ip'. For the mo- 
ment we ignore the magnetic field and assume a real GL function so that 
the GL equation reads £ 2 V 2 V> — tp + V' 3 /'0 2 <) = 0. Since tp' is proportional to 
1/6, we keep its linear terms only, £ 2 V 2, 0' + 2ip' = 0. This equation has the 
exponential solution ip'{x) = —ip oc c~^ x ^£ i /(\/2b). Since the GL coherence 
length £ is much smaller than the extrapolation length b, the exponential is 
small compared to the constant term. Assuming that the GL function pro- 
vides us with the order of magnitude estimate of the gap function, we find 
that 5c = 5A(0)/A w 5ip(0)/ip = -(£/V2)6(l/b). After substitution of this 
estimate into Eq. (1.14) we obtain 



1 A(0) 1 1 5N^ 
1 + -^ A e(0)N V N • 



(1.15) 



We can assume, that the induced density of states per area is linearly 
proportional to the applied electric field, 



6NW =Eg 



(1.16) 



1 Surface superconductivity controlled by electric field 



7 



so that the modified inverse length is of the form 

where U s has a dimension of a potential. According to (1.15) this effective 
potential is given by 

wr- v m^-k- (L18) 

We have introduced a dimensionless parameter 

which captures the effects of the gap profile near the surface. According to 
de Gennes' estimate [5] , the surface ratio r\ is of the order of unity. 

Let us note, that a boundary condition similar to the de Gennes boundary 
condition described above can be derived also from the minimum free energy 
principle [12]. 



1.4 Link to the limit of large screening length 

To draw a link between the de Gennes-type formula (1.18) and the field- 
affected GL boundary condition (1.8) obtained in the limit of large Thomas- 
Fermi screening lengths, we evaluate the coefficients of the de Gennes formula 
in the jelly model. As above we assume that in the absence of the electric field 
the extrapolation length diverges, l/bo = 0. 

For zero electric field the density of states is step-like, N(x) — Nq for x > 
and N(x) = elsewhere. Now we include the electric field. It is exponentially 
screened due to the induced electron density given by (1.1). We employ the 
local density approximation and assume that the local density of states is a 
function of the local density of electrons 

N(x) = N + ^5n(x). (1.20) 
on 

This approximation yields a simple change of the density of states per area 

oo 

SN^ = ^ J dx6n{x). (1.21) 
o 

The induced density of electrons per area is given by the Gauss law 

oo 

eE = -e J dxSn{x). (1.22) 
o 
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From relations (1.16), (1.21) and (1.22) we find the coefficient of the density 
of states 

g= -^e (L23) 
To be able to compare the BCS formula with the de Gennes-type one, we 
have to express both expressions in terms of the same parameters. We thus 
convert the parameters of de Gennes-type formula to their phenomenological 
counterparts. 

The GL coherence length £ = h/^—2m*a depends on the temperature 
according to £ = £(0)/a/1 — T/T c , where £(0) is the 'zero' temperature value. 
From a = a'(T — T c ) one finds 

— - (1.24) 



em n 

Finally, we have to express the BCS interaction potential V in terms of 
the critical temperature. In the BCS critical temperature 

T c = 0.85 6>dc _1 ^ (1.25) 

we assume that the Debye temperature 6>d and the BCS interaction V are 
independent of the electron density. This corresponds to approximations we 
have tacitly used above ignoring the electric field effect on the phonon spec- 
trum. The density dependence of the critical temperature follows thus from 
the density dependence of the density of states 

dT c i 1 8N Q 1 dN a 

— =O.850 D c ^^—=1^-^. (1.26) 

We will use this relation to express density derivatives of the density of states 
in terms of the density derivative of the critical temperature. 

Now we can rewrite the effective potential in terms of phenomenological 
parameters. Using £(0) from (1.24) in equation (1.18) we find 



2m*a 

u s ■' IT • \ \ 

Next we substitute g from (1.23) 



- = -V^^T c -^-g. (1.27) 



1 2m* a' 1 8N e 

Trr"—^^ (L28) 

The group of terms around 8N /dn can be substituted with the help of (1.26) 
so that we obtain 

1 2m*a'dT c e 

wr^^r^-e (L29) 

The boundary condition we have obtained now from the de Gennes-type 
formula 
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dip E N 2m*a'dT c eE 

differs from the large screening length limit (1.8) by the factor r\. Of course, 
a heuristic derivation of the field-effect from the GL equation cannot cover 
the factor r\ which depends on the gap profile on a scale smaller than the GL 
coherence length. 

In summary, the electric field applied to the surface of the superconductor 
modifies the GL wave function near the surface. This effect is conveniently 
described by the GL theory, where the GL equation remains unaffected by 
the field and the entire electric field effect is covered by a modified boundary 
condition. In the next section we discuss an experiment which can be used to 
measure the predicted field effect on the GL boundary condition. 



1.5 Electric field effect on surface superconductivity 

We will investigate now the magneto-capacitance for magnetic fields near the 
surface critical field B c ^. We focus on this region since we expect that the bias 
voltage affecting only the surface has a relatively large effect on the surface 
superconductivity. In this section we show how the electric field affects the 
nucleation of superconductivity [20] . 

1.5.1 Nucleation of surface superconductivity 

At the surface critical field B c ^ the superconductivity nucleates in the surface 
region. At the nucleation point the GL wave function ip is infinitely small, 
therefore we can work with the linearized GL equation, omitting the cubic 
term in the equation (1.2) 

— ^(-ifiV - e*A) 2 V> + cm/; = 0. (1.31) 

The solution is restricted by the boundary condition (1.9). 

The electrode is a superconductor which fills the half space x > 0. We 
assume a homogeneous applied magnetic field B a = (0,0, B a ). Since an 'in- 
finitely' large electrode has translation invariance along the y direction, we 
use the Landau gauge of the form 

A = (0,B B a;,0). (1.32) 

Nucleation is possible if (1.31) has a nonzero solution, i.e. if the pa- 
rameter — a becomes equal to an eigenvalue e of the kinetic energy given 
by — (—iftV — e*A) ip = sip. Since a changes with temperature, a = 
a'(T — T c ), the eigenvalue e of the kinetic energy determines the nucleation 
temperature T* according to T* — T c = —e/a'. To avoid dual notation for 
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the same quantity, we will treat the equation as an eigenvalue for a. Since a 
is negative, the nucleation temperature T* is always below the critical tem- 
perature T c in the absence of the magnetic field. Note that we are looking for 
maximal a. 

Assuming the translation invariance along the y and z axes we can write 
the wave function as 

1>(x,y,z) = ?p(x) e lky e iqz . (1.33) 
Using (1.33) in the GL equation (1.31) we get a one-dimensional equation 

s(-(5)' + ( t -X , )' + «')* + "*- a < 1JM > 

Any non-zero value of q results in the kinetic energy q 2 h 2 /2m which lowers 
the value of a reducing the nucleation temperature. The nucleation happens 
at the first possible occasion, i.e., at the highest allowed temperature. We thus 
take q = 0. Similarly, we have to find the wave vector k from the requirement 
of the highest nucleation temperature. 



1.5.2 Solution in dimensionless notation 

It is advantageous to express the x-coordinate with the help of the dimension- 
less coordinate r 

x = rl + 2l 2 k, (1.35) 

with the magnetic length 

I 2 = — ^_ (1.36) 

and the momentum to = —2kl. The wave function is then proportional to the 
parabolic cylinder function of Whittaker [18] 

^(x)=ND (j+to), (1-37) 

which solves the differential equation (1.31) in the dimensionless notation 

The dimensionless boundary condition is 



A,(r ) 



T = -2kl 



fV"+5' (1-39) 



where the prime denotes a derivative with respect to r . The parameter 

1 am* 

v = -~ TTT (1-40) 

2 e*hB a v ' 
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plays the role of an eigenenergy. 

The boundary condition (1.39) links v and To- Since we are looking for 
the minimal v, we take the solution of (1.39) as a function v{t ). Besides the 
obvious numerical search we can give directly a nonlinear equation for this 
desired minimum given by v'{tq) — 0. For this purpose we differentiate (1.39) 
with respect to To arriving 5 at [20,21] 



D p+1 (t ) 




Ato) . =- 2 ^ + i )(1+ ||) V ^ ' 2 J ^ ' 5 V bV 1 2' 

(1-41) 

where v is the minimal value of v. 

Solving equation (1.41) we find the minimal v to each given to, i.e., 0[to]. 
Since to = —2kl, we find in this way the eigenenergy as a function of momen- 
tum k. The GL wave function (1.37) is now specified except for its amplitude 
TV. We discuss this amplitude below. In figure 1.1 we see how the shape of 
the GL wave function evolves with inverse extrapolation length 1/6, i.e. how 
it depends on the external bias. For positive electric fields attracting charge 
carriers to the surface, the superconducting density is pushed from the surface 
into the bulk while for oppositely directed electric fields the superconductivity 
is even more squeezed near the surface. 

The lowest eigenvalue v corresponds to the highest attainable critical mag- 
netic field 

-m*a B C 2 



*c3 = ^t^tt - 7TFTTT' (1-42) 



where B c i is the upper critical field. In figure 1.2 we present the result for the 
surface critical field (1.42) versus the external bias. Without external bias the 
known GL solution [22] with B c3 /B c2 = 1.69461 is reproduced. We see that 
the external bias can enhance or decrease the surface critical value. 

According to the GL wave function (1.33), the current flows only in the y 
direction. Its net value given by the x integral of the current density equals 
the k derivative of the eigenenergy. Thus, for the minimal v the net current 



5 With the minimum condition v' — at To, the derivative of (1.39) yields 



D- 



D'- 2 



The left hand side follows from (1.38) as D"/D„ = ?f - v - \. The right hand 
side is given by (1.39). The result is 



t = ±2W 1 + 



with the negative root being the physical one. Finally, substituting this to and a 
general relation D'^ — ToDp/2 — D D +i into (1.39) we obtain (1.41). 
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Fig. 1.1. The GL wave function of the condensate versus the boundary condition 
(1.9). The normal GL wave function without external bias is marked as thick line. 

is zero. Current distributions for three different boundary conditions with 
extrapolation lengths £/b = 0, +1 and -1 are plotted in the fig. (1.3) as bold, 
dashed and dash dot lines. Since there is no net current circulating around the 
sample, as is seen in the inset of the fig. (1.3), the magnetic field is reduced 
only in the region of nucleation. Each surface thus acts independently. 

1.5.3 Surface energy 

With the help of the GL wave function we can now calculate the surface 
energy 



Here we neglect the free energy of the magnetic field, but the term ~ \ip\ 4 
which determines the amplitude of the GL wave function cannot be omitted 
even near the nucleation line. 

Since the applied field changes the shape of the GL wave function only 
weakly while its amplitude changes rapidly near the critical point, we fix the 
shape to be the nucleation function (1.37) and use the amplitude TV as a 
variational parameter. Briefly we substitute tp = AfDp into (1.43) which gives 



oo 



a = 




[B a -B(x)} 2 | |(iftV + e*AM 
2/zq 2m* 



(1.43) 



o 



a = 




(1.44) 



where a — — (Tie*B/2m*)(2D(E,B) + l) is the maximal eigenvalue corre- 
sponding to the nucleation temperature under given magnetic and electric 
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Fig. 1.2. The surface critical field B C 3 versus extrapolation length. The exact solu- 
tion (1.42) is the solid line and its slope at £/b = is shown by the dotted line. The 
inset shows the ratio of integrals (1.45) and (1.46). 



field, while a — a'(T — T c ) is the GL parameter given by the actual sam- 
ple temperature. The quadratic and quartic terms in (1.44) are weighted by 
dimensionless integrals 



I 2 = I drDl(r), 

OO 

h = J cItDUt). 



(1.45) 
(1.46) 



In the normal state, a > a and the minimum of the surface energy is Af 2 = 
giving a = 0. In the superconducting state a < a and the minimum of a 
(1.44) is at TV 2 = —(a — a)/ 2 /(/3^4) giving the surface energy 



h 2(3 



(1.47) 



At the critical point a = a and the surface energy vanishes. The surface 
energy and its first derivatives with respect to the electric and magnetic field 
are continuous at the critical point. The second derivatives are discontinuous. 
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2 4 „ 6 8 
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Fig. 1.3. Distribution of the current density for boundary conditions with extrap- 
olation lengths £/6 = 0, +1 and -1 are plotted in bold, dashed and dash dot lines. 
The inset shows the corresponding reduction of magnetic field. 



In the next chapter we show how this discontinuity appears in the magneto- 
capacitance. 



1.6 Magneto-capacitance 

We assume a capacitor in which the first electrode is a superconducting and 
the second electrode is a normal metal. Our aim is to evaluate the contribution 
of the surface superconductivity to the capacitance. 

The capacitance of the capacitor with one superconducting electrode reads 

where S is the area of the capacitor, C n is the capacitance when both elec- 
trodes are normal, e is the permittivity of the ionic background in the su- 
perconductor. This follows from the inverse capacitance given by the sec- 
ond derivative of the total energy W with respect to the charge Q, 1/C = 
d 2 W/dQ 2 . The charge Q is linked to the electric field E at the surface of the 
superconductor via the Gauss law eE = Q / S. Since energies of normal and 
superconducting states differ by the surface energy W s = W n + Sa, one arrives 
at equation (1.48) for the difference in capacitances. 
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1.6.1 Discontinuity in magneto-capacitance 

Now we can evaluate the jump of the capacitance, which appears as the mag- 
netic field B exceeds the critical value B c3 . Since as — > a for B — ► £? c3 , the 
discontinuity of the inverse capacitance equals 



1 1 n 2 e* 2 B 2 3 l I 2 ( d, 



- 



C s C n e 2 m* 2 S(3I 4 \dE 



(1.49) 



where we have used daE/dE = — (he*B/m*){dv/dE). 

To evaluate the slope dv/dE we recall the numerical result shown Fig. 1.2. 
The tangential dotted line yields B c3 /B c2 = l/(2z>+l) = 1.69- 1.69^/6 which 
follows also from an explicit variational calculation, see Eq. 20 in [20]. Since 
d(l/b)/dE = l/U s , for 1/6 -> we find {dv/dEf = 0.087 (?/U*. From 
relation (1.29) thus follows 

0-087-4— -CW hnr^ • (i-so) 



dE ) ' h 4 (e*) 2 V 9 Inn 

We have used the relation [3 = a'T c /n which holds for Gor'kov values of GL 
parameters. Substituting (1.50) into (1.49) we arrive at 

JL - J_ - -0 348 ^^£W (1 51) 

The GL coherence length £ = hj\j2m*a relates to the surface critical field, 
see (1.42)). For 1/6 — > it yields £ 2 = 1.69fi/(e*B C 3). Moreover, according to 
(1.36) we can express B c3 via the magnetic length, therefore 

We have used Jf/Zi = 2.42, which is the value at 1/6 = 0. 

Since the GL parameter (3 can be fitted from experimental results, rela- 
tion (1.52) allows one to establish r/ (d\nT c /d\nn). This material parameter 
describes the change of the critical temperature with the electron density. 



1.6.2 Estimates of magnitude 



For an estimate we assume some typical numbers. The most sensitive measure- 
ments of capacitance performed in the C ~ ^iF range are capable to monitor 
changes SC/C <~ 10~ 6 with error bars at SC/C ~ 10~ 7 . From the capacitance 
C = edS/L one sees that a 1000 A-thick dielectric layer with = 10 3 eo has an 
optimal area of S = 10 mm 2 which is about the usual size of such samples [10]. 

To estimate (3 we use Gor'kov's relation (3 = 6n 2 k^T 2 / (7((3)Epn). For 
Niobium T c = 9.5 K and n — 2.2 x 10 28 /m 3 . The free electron model used by 
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Gor'kov than gives the Fermi energy Ep = 4.6 1CP 19 J. The corresponding 
GL parameters then is (3 = 1.2 10 -53 Jm 3 . The logarithmic derivative is 
estimated in [14] as dh\T c /d\nn = 0.74. Since r] is not known, we take r] = 1 
according to the simple theory. Finally we need the third critical magnetic 
field B C 3 to estimate the magnetic length I. From B c3 = 1.69 B C 2 and the 
experimental value B C 2 = 0.35 T [6] one finds £? C 3 = 0.59 T, which yields 
I = 325 A. With these values from equation (1.52)we obtain the discontinuity 
1/C S — 1/C n = —2.8 10~ 4 F _1 . Since the capacitance was estimated to be 
<~ 10~ 6 F, the corresponding relative change C s /C n — 1 <~ 3 x 10~ 10 is too 
small to be observed. A slightly more optimistic estimation can be found in 
Ref. [21]. 

In high-T c materials the GL parameter (3 is by three to four orders of mag- 
nitude larger than in conventional metals due to larger T c and lower density 
of holes giving also lower Fermi energy. Moreover, the logarithmic derivative 
of the critical temperature is about <91nT c /<91nn = —4.82 as estimated in 
Ref. [14]. Finally, higher critical surface magnetic field B c3 allows to reduce 
the magnetic length to values limited rather by experimental facility. For a 
typical field of 10 T the magnetic length is 79 A. These factors together provide 
an enhancement to values 1/C S — 1/C n = —52 F _1 or C s /C n - l~5x 10~ 5 
which is an experimentally accessible discontinuity. 

1.7 Summary 

We have shown that the electric field applied to the surface of the supercon- 
ductor modifies the boundary condition of the GL wave function. Since the 
surface superconductivity is sensitive to this boundary condition, we have dis- 
cussed the influence of the electric field. From the surface energy we predict 
that a planar capacitor with one normal electrode and the other electrode 
to be superconducting reveals a discontinuity of the capacitance at the third 
critical field B c3 . This discontinuity is too small for capacitors from conven- 
tional superconductors but it is large enough to be observed in capacitors 
with ferroelectric dielectric layers of a width of 1000 A and non-conventional 
superconductor electrodes. 
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